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An automatic frequency control (AFC) loop is introduced and analyzed in detail The 
new scheme is a generalization of the well-known Cross-Product AFC loop that uses run- 
ning overlapping discrete Fourier transforms (DFTs) to create a discriminator curve . 
Linear analysis is included and supported with computer simulations. The algorithm is 
tested in a low carrier to noise ratio (CNR) dynamic environment , and the probability of 
loss of lock is estimated via computer simulations. The algorithm discussed is a subopti- 
mum tracking scheme with a larger frequency error variance compared to an optimum 
strategy , but offers simplicity of implementation and a very low operating threshold CNR. 
This technique can be applied during the carrier acquisition and re-acquisition process in 
the Advanced Receiver. 


I. Introduction 

Automatic frequency tracking is a subject of great impor- 
tance in the fields of communications, control, and signal 
processing. For instance, it may provide initial rapid acquisi- 
tion in coherent receivers, may supply information on the 
velocity of rapidly moving targets, and can be used in con- 
junction with channel doppler extractors, just to name a few 
applications. 

This article presents a study of the fundamental problem of 
estimating the frequency of a sinusoidal wave embedded in 
noise, subject to severe dynamics. It is well known [1] , [2] 
that the maximum likelihood estimator (MLE) of a constant 
(but unknown) frequency sine wave is equivalent in practice 
to a discrete Fourier transform (DFT) operation on the re- 
ceived discrete data, the MLE of frequency being the center 
frequency of the bin filter in the DFT with maximum output 


power. For a stable frequency, arbitrarily good resolution is 
obtained by increasing the size of the DFT or, equivalently, 
by extending the observation interval. However, when the fre- 
quency of the incoming wave is time varying in a random man- 
ner, the implementation of a maximum likelihood approach 
may be computationally demanding, as discussed in [3] , [4] . 

The object of this article is to propose and analyze a sub- 
optimum frequency control loop based on short overlapping 
DFTs and suited to track the rapidly varying frequency of a 
sine wave. For brevity, the new scheme will be referred to 
here as ODAFC (Overlapping Discrete Fourier transform-based 
Automatic Frequency Control). This scheme can be thought 
of as a generalization of the so-called Quadri-Correlator [5] . 
The new scheme integrates ideas originating in the fields of 
classical spectral estimation and digital phase-locked loop 
theory. It will be shown that the ODAFC is capable of track- 
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ing the frequency of a carrier signal in a low carrier-to-noise 
(CNR) and high dynamic environment where a Phase-Locked 
Loop (PLL) is inoperative. Such a situation is foreseeable dur- 
ing planetary encounters. In these cases ODAFC will allow 
fast carrier re-acquisition and help the resumption of tracking 
itself. This is achieved with remarkably modest implementa- 
tion complexity, since the algorithm may be implemented in 
software and does not require any modification to the Ad- 
vanced Receiver (ARX) hardware. 


are subtracted every sampling interval (a weighted combina- 
tion of the filter bin outputs was not selected to facilitate the 
implementation). The subtraction of the powers in the bin 
filters around zero frequency to create a discriminator curve 
as opposed to selecting the bin filter with maximum power 
was originally discussed in [7] . The idea of employing run- 
ning overlapping DFTs has a connection with the use of the 
FFT for the estimation of power spectra based on averaged 
periodograms [8] . 


First, as an indicator of the quality of the closed loop esti- 
mator, the variance of the loop is computed for a fixed fre- 
quency sine wave in the presence of additive white Gaussian 
noise (AWGN). Later, the loop is tested with severe dynamics 
and noise together, and the probability of loss of lock is esti- 
mated. Both analytical answers (whenever possible) and com- 
puter simulation results are included for greater confidence. 

II. Description of the ODAFC Loop 

A. System Block Diagram 

A block diagram of the ODAFC loop is shown in Fig. 1(a). 
An equivalent baseband model is shown in Fig. 1(b). The dis- 
crete time in-phase and quadrature mixer outputs admit a 
representation of the form 


III. Linear Tracking of the ODAFC Loop 
in the Presence of Noise 

A. Exact Computation of the Variance of the 
Frequency Error 

This section focuses on the tracking performance of the 
ODAFC loop when the carrier-to-noise ratio is sufficiently 
high to justify a linear analysis. A sine wave with constant fre- 
quency is assumed for the remainder of the section. The goal 
is to compute analytically an expression that predicts the 
variance of the closed loop estimator. The results will be con- 
firmed by computer simulations. 

The closed loop estimator uses the received complex sample 


/ = A cos (4> n ) + n In 

(la) 

n n 

Q n = A sin (4> n ) + n Qn 

and computes 


(lb) 



The amplitude A is related to the carrier power via A 2 -P c , the 
instantaneous phase error is given by <p n (rad), and the noise 
sequences {n In }, } are independent Gaussian random 

variables with zero mean and variance o 2 =N 0 /2T S , where T 
is the loop update time (henceforth called the sampling inter- 
val). The carrier-to-noise ratio is defined as CNR = PJN 0 . 
Every update time, N s complex samples are Fourier trans- 
formed with N s zeroes appended. At the tick of the clock, the 
data are shifted when another complex sample arrives; a new 
discrete Fourier transform is computed (of size 2N S ), an error 
control signal is created and the loop filter is updated. This 
process is repeated every sampling interval T s , each time incor- 
porating just one more complex sample. A timing diagram of 
the underlying process is illustrated in Fig. 2. 
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The discriminator output at time n-k is given by 
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B. On the Discriminator Characteristic 

It is well known that the discrete Fourier transform can be 
seen as a bank of bandpass filters tuned at multiples of half the 
Nyquist rate [6] . In order to create a control signal propor- which upon ex P ansion results in 
tional to the frequency error (discriminator characteristic), 

the powers in two adjacent filter bins around zero frequency ? k - ^ + fV q k 
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where 


center frequency is moving, since larger DFTs have a smaller 
linear operating range. 


S* = 2 (f) E E A 2 & 

\ sl m =k-(N s -l) n=k-(N s -l) 

X sin ^ (m - «)j 

= 2 ljrj £ £ 

V s/ m =k-(N„- 1) n=k-(N c -l) 


The term N eq k will be referred to as the additive equivalent 
noise. Unfortunately, the typical white noise assumption for 
the equivalent noise is not justified here, since there exists a 
significant amount of correlation between the noise samples. 
It is shown in Appendix A that the correlation function of the 
noise sequence {Af eq ik } is given by 


R (8) = E (7V eq ,fc ^eq ,k+fi) 

iV eq 


+ A n i„ s*n (<P m )~An [m sin (0„) 

- A n Qn C0S + H ln n Qm " U Qn n im j 

X sin|^- (m - «)j ( 7 ) 


for a fixed frequency error; S k ( •) denotes the discriminator 
characteristic for the ODAFC loop, and it is given more explic- 
itly in Eq. (8). 



The slope at the origin is 


S^ k T s ) 


Aw fc 3 r° 
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otherwise 




R n (£) = 0 

Jv eq 


The discriminator characteristic is illustrated in Fig. 3 for a 
unit amplitude signal ( A = 1) for N s = 2,4 corresponding to 
DFTs of sizes 4 and 8, respectively; notice that a larger slope 
for the discriminator at the origin is obtained with larger 
DFTs. For a fixed frequency sine wave in the presence of 
AWGN, the larger the size of the DFT, the smaller the variance 
of the estimator. However, in a dynamic medium, the maxi- 
mum permissible size for the DFT is limited by how fast the 


With the aid of Fig. 1 , it is a straightforward exercise to show 
that in operational notation the Z transform of the normalized 
radian frequency error is given by 

Afi(z) T = H{z) O (z) T - ^ N (z) (1 1) 

* S'JO) 
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where 


S’(0)F(z)NCO(z)k N 

H(z) £ — (12) 

1 +S^(0)F(z)NCO(z)k N 


represents the closed loop transfer function. The loop filter 
F(z) considered here has the form 




1 -Z' 1 


Unfortunately, the previous integral does not reduce to a 
simple compact expression as in a standard phase-locked loop 
[5] , and a numerical integration must be used. For illustration 
purposes, Fig. 4 shows numerical results, based on Eq. (15), 
for N s = 4 and for a particular mechanization using T s = 2 X 
10" 3 s with nominal bandwidth ^(Hz) as a parameter. Com- 
puter simulation results are also included in the same graph to 
confirm the analysis. Notice that for high CNR, the analysis 
and the simulations are in excellent agreement. 

B. A Useful Approximation for the Variance of the 
Frequency Error 


1 -z 






where 


k 


l 


r + 1 



r 


A very useful approximation to the variance of the fre- 
quency error can be obtained for low to moderate CNR when 
the bandwidth of the loop is much smaller than the sampling 
rate. The approximation is obtained by expanding the cosinu- 
soidal terms of the noise spectral density of Eq. (16) in a 
Taylor series and keeping only the significant terms of the 
expansion. The result is 



an L T s~^L T s) 3 


(17) 


r = 4 £ 2 


where 


£ = damping ratio 
B a = nominal bandwidth (Hz) 

and NCO{z) is the mathematical model of the numerically 
controlled oscillator and may include computational delays 
inherent in a digital implementation. Of particular interest is 
the NCO discussed in [9] , which conceptually represents the 
cascading of an integrator using the trapezoidal rule and two 
extra delays. 

r(z + 1) 

NCO(z) = — (14) 

2 z 2 (z - 1) 

Since the noise sequence {A£ is stationary, the steady state 
variance of the frequency error (Hz) is given by [10] , [11] as 



= 2tt( 2B l ) 

B l = one-sided loop noise bandwidth 

The approximation is very tight for small bandwidths 
(. B l ^ 10 Hz) and for low CNRs (^ 30 dB-Hz), and it loosens 
as these two parameters increase, until it becomes useless for 
CNR y 50 dB-Hz, as shown in Fig. 4. 
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where the equivalent noise spectral density has the form 


IV. Steady State Frequency Errors 

The steady state frequency error due to acceleration in fre- 
quency can be calculated from Eq. (1 1) in the absence of noise 
using the final value theorem [12] . The input frequency in the 
Z domain is given by 


and the steady state error is 

A4(Hz) = 

_ 'o('+l) 2 1 

16 B?r 2n 
where the units of the jerk/ 0 are in rad/s 3 . 

V. The ODAFC Loop in a Dynamic 
Environment 

One of the main reasons for this study was the need to 
operate the frequency estimator under severe dynamic condi- 
tions at low CNR levels. It is natural to select the probability 
of loss of frequency lock as a measure of quality. Once cast 
into this form, computer simulations provide the most direct 
way to evaluate performance. 

A. Dynamic Trajectory 

A dynamic trajectory of duration 8 s was developed to 
evaluate loop performance. Initially, it is assumed that the 
phase, frequency, and frequency rate of the received signal 
are perfectly known, i.e., that the loop is in the tracking mode 
of operation. The transmitting vehicle then experiences a sud- 
den maneuver resulting in an acceleration of 3 s duration and 
equivalent to a doppler rate of -1287 Hz/s. For the next 0.5 s, 
a frequency acceleration of 5150 Hz/s 2 occurs, followed by a 
2 s duration doppler rate of 1287 Hz/s. A frequency accelera- 
tion of -5150 Hz/s 2 proceeds for the next 0.5 s, and the 
maneuver culminates with a 2-s-duration doppler rate of 
-1287 Hz/s. 

B. Probability of Loss of Lock for the ODAFC Loop 

The probability of losing frequency lock was estimated for 
a loop that is initially in-lock and then experiences the phase 
trajectory previously described. The simulations were carried 
out for various sizes of DFTs, but N s = 4 was found to have the 
lowest possible threshold with this technique. Some simulation 
results are presented in Figs. 5 and 6 for A^ = 4 and 8, respec- 
tively, as a function of carrier to noise ratio (dB-Hz) with 
nominal loop bandwidth as a parameter. For a given CNR 
and nominal bandwidth, 250 different simulation runs were 
employed. Loss of lock was declared when the instantaneous 
frequency error exceeded the one-sided Nyquist bandwidth 

1/2 T s . 


In order to put these results in proper perspective, from [4] 
a type-3 digital phase-locked loop (DPLL) with the same phase 
trajectory and an optimized bandwidth, has a probability of 
0.1 of losing lock at approximately 26 dB-Hz; i.e., the thresh- 
old for a DPLL is about 3 dB higher than that of the ODAFC 
with N s = 4. From the same reference, the threshold of a typi- 
cal Cross-Product AFC loop is 2.3 dB-Hz higher than that 
of the ODAFC. 

Estimation of the probability of loss of lock is only a par- 
tial performance measure. If in-lock, a phase-locked loop 
has typically a better frequency estimation ability than the 
ODAFC loop. In Figs. 7 and 8, the rms value for the ODAFC 
loop is estimated for the aforementioned trajectory. Notice 
that when N s = 8, smaller tracking errors are achieved, but the 
threshold is higher. The same effect is noticed for larger values 
of N s . This scheme then suffers from the usual trade-off be- 
tween noise suppression and dynamic tracking: the effect of 
the noise increases as N s decreases, while smaller dynamic 
errors are obtained as A^ decreases. 

VI. Conclusions 

A new automatic frequency control loop was introduced 
and discussed. The algorithm (ODAFC) is based on running 
overlapping DFTs, and can be viewed as a generalization of the 
well-known Cross-Product AFC (CPAFC) loop [7]. In fact, 
the CPAFC is a special case of the ODAFC loop when A^ = 2. 

A detailed noise analysis was presented of the new algo- 
rithm for any desired DFT size, and the analysis supported 
with computer simulations. In addition, the algorithm was 
tested in a noisy dynamic medium, and the probability was 
estimated of loss of lock for various configurations for a spe- 
cific phase trajectory. It was shown in the text that for this 
particular trajectory, a value N s = 4 minimizes the operating 
threshold for this algorithm. In fact, the ODAFC with A^ = 4 
has a threshold 3 dB-Hz lower than that of a DPLL and com- 
parable to that of a maximum likelihood estimator [4] , but 
with a larger error variance. The implementation complexity 
for the ODAFC is substantially simpler than for an optimum 
strategy, not even requiring fast Fourier transforms to com- 
pute the DFT when A^ is small. 

In summary, the algorithm discussed in this article is a 
suboptimum tracking scheme with increased variance for the 
frequency error compared to an optimum strategy, but offers 
simplicity of implementation and a very low threshold. 
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Fig. 1. Two loops: (a) an automatic frequency control loop using 
overlapping DFTs and (b) an equivalent linearized baseband model 



FREQUENCY ERROR, Hz 


Fig. 3. Normalized discriminator curves for ODAFC loop (only 
positive frequency errors are shown) 
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Fig. 2. Overlapping of the received signals to create a 
discriminator characteristic 



CNR, dB-Hz 

Fig. 4. Tracking performance of the ODAFC loop ( Af # = 4) as 
a function of CNR with nominal bandwidth as a parameter 
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g. 5. Probability of losing frequency lock in a dynamic medium as Fig. 7. Tracking performance of the ODAFC loop in a dynamic 

function of CNR with nominal bandwidth as a parameter (N s = 4) medium (N s = 4) 
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Fig. 6. Probabilityof losing frequency lock in a dynamic medium as Fig. 8. Tracking performance of the ODAFC loop in a dynamic 

a function of CNR with nominal bandwidth as a parameter (N $ = 8) medium (N s = 8) 
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Appendix A 

Correlation Function of the Equivalent Additive Noise for the ODAFC Loop 


In this appendix, the autocorrelation function of the Due to symmetry, the previous correlation reduces to 
equivalent additive noise for the ODAFC loop is computed for 
any lag. The starting point is Eq. (7). The autocorrelation is 
given by 


R m = E(AT AT ) 

eq 
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For linear analysis it is useful to make the approximation 
cos (<fi x - <t> y ) = 1 in Eq. (A-l); therefore, the correlation func- 
tion becomes 
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or, using short notation: 
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